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We compare the Radon transform in its standard and symplectic formulations and argue that 
the inversion of the latter can be performed more efficiently. 
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1. Introduction 

The Radon transform (Q) is a key mathematical tool in tomography. Its inverse 
enables one to reconstruct a function if some of its integrals are known. The whole 
subject has been recently revived by quantum mechanical applications. The pos- 
sibility of reconstructing the tomographic map of the Wigner quasidistribution 
function SQl associated with a given quantum state (0; Q; B) has motivated 
experiments jlHl) triggered novel proposals (H ) and boosted innovative the- 
oretical techniques (jXi). Applications are widespread and diverse. The entire field, 
driven by a blizzard of technical advances, is attracting increasing attention and is 
growing at a lively pace. Good reviews on the subject can be found in (Idl). 

The Radon transform was originally introduced as an integral transform defined 
over submanifolds of IR", that may be viewed as a "configuration space." However, 
if n is even, one may think of IR" as a phase space and consider the integrals over 
its Lagrangian submanifolds. One may then associate the tomographic map with 
the symplectic transform on the phase space (Q)- I^i this context, motion is instru- 
mental for the identification of the phase space and its Lagrangian variables: the 
Hamilton equations do not appear in the definition of the Radon transform and this 
interpretation differs from the original one. Nevertheless, the approach is prolific 
and enables one to identify different types of tomograms (0), extend tomography 
to curved surfaces (fTB ) and consider more general problems and applications. This 
is in line with previous hystorical developments, by Radon himself 0), John (0), 
Helgason (fi^ ) and Strichartz IH), and paves the way towards so far unhearted 
quantum mechanical applications. 

In this article we shall compare the standard Radon approach with that based 
on the afore- mentioned symplectic identification. We shall argue that, although 
mathematically equivalent, they may differ in practice. In particular, the inversion 
may be far from trivial and may turn out to be simpler in the symplectic framework. 
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2. Symplectic tomography 

Let us focus on the 2-dimensional case for the sake of concreteness. The Radon 
transform, in its original formulation, solves the following problem: reconstruct a 
function of two variables, say /(p, g), if its integrals over arbitrary lines are given. 
The Radon transform (or homodyne tomogram) reads 

f{^,X)= [ f{q,p)S{X -qcos^-psm^)dqdp, (1) 

where 5 is the Dirac function, G [0, 27r), to = (cos "i?, sin??) G S (the unit sphere in 
ID) and X £ R. In order to obtain a symplectic formulation, a central observation 
is the following: it is possible to express the Radon transform in affine language 
(the so-called tomographic map) (Q, [13) and write 

f^'" (/i, u,X)= [ f{q, p)6{X -fiq- vp) dqdp, (2) 

where p,,v,X S IR, + v"^ > 0. We have named "M^" the tomographic map ([2]) 
after Man'ko and Marmo, who gave seminal contribution towards its significance 



\M,\M,m,\2i). Clearly 

/«(??,X) = /^^'(cos??,sin??,X). (3) 

Consider now a particle moving on the line g E IR and a function f{q,p) on its 
phase space {q,p) G IR^. Since 

(''^)(j?)(:)--+--(--)(roO(:)' 

the argument in the Dirac delta function in Eq. ([2]) may be considered either 
as a Euclidean product or as a symplectic product. The two interpretations are 
completely equivalent and one can equivalently solve the inversion problem by 
using the Euclidean or symplectic Fourier transform. We shall see in the next 
section that the two procedures can vastly differ in complexity. 
Note that the Radon transform is defined in an equivalent way by 

f\^,X)= [ f{scos^-Xsm^,ssm^ + Xcos^)ds. (5) 
Jr 

The inversion formula, as given by Radon, amounts to consider first the average 
value of /" on all lines tangent to the circle of center {q,p) and radius r, namely. 



1 f 

— / f\'&, q cos 'd + p sin + r)d'i!} (6) 
27r Jo 

and then compute 



where F'^^ (r) denotes the derivative with respect to r. The Radon transform maps 
a (suitable) function on the plane into a function on the cylinder. Some conditions 
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that guarantee the invertibihty and continuity of the map were studied by Radon 
himself Q), John (B[), Helgason (H) and Strichartz HI). 

On the other hand, the inverse transform of ([2|) reads (f^ . [H) 



3. The inverse transform: an explicit example 

We now compare the inversions ©-([T]) and ^ by looking at a very simple example: 
the ground state of a one dimensional quantum harmonic oscillator, 



vr 

Its M^-transform reads 

„2 



f{q,p) = -e-'''(i'+P'-l (9) 



'2 /» 

f^'\n,u,X) = — e-"'('''+^^')(5(X -nq- up) dqdp 



a2 



a2 



dp 



= , " e-"-^, (10) 
^7r(/x2 + 1/2) 

which is a Gaussian with respect to X, but has a nontrivial dependence on \i and 
I/. On the other hand, by making use of dS]), one gets the Radon transform 

f\^,X) = f^\cos^,sm^,X) = ^e^"'^', (11) 



which is simply a Gaussian, independent of the angle 'd, due to symmetry. 

Let us tackle the inversion problem. We start from the inverse transform, 
which is easily solved in a few lines: 

fiQ^P) = 77^ I f^\^^,iy,X)e'(''->^'^-''PUxdfidu 



1 



,,2 , ,,2 



IR2 



(2^)2 

(12) 



Let us now endeavour to invert the Radon transform. It would be tempting to 
leave this as an exercise for the reader, but we will sketch the main steps of the 
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derivation. From ([6]) we get 

F' (r) = — ^/ (gcosi9 + psini9 + r)e-"'('?=°'^''+P'^i'^''+'')'di? (13) 
vrVvr Jo 

and thus f{q,p) = Yvoieio fe{q,p), where 

Already in the very simple case of a Gaussian function with a Gaussian Radon 
transform the above inversion formula is not easy to manage. First introduce a step 
function, 9{r) = 1 only if r > 0, and change the period of the angle integration 

f^{q,p) = -^ [ dr r d'd ^^'' ~^\ r-qcos^+ -psin^^je-"'^"-'?^"'^^-^'^"^'^)'. 
vt^VtF r 

(15) 

Then change variables z = r — qcos-d — psm.'d and t = tan(T9/2) 

„, , 2a^ ('(t^{z-q-e) + 2pt+{z + q-e 



T^'^V^ Jr^ V 1 + 

dtdz. (16) 



ze 



t'^{z-q) + 2pt + (z + : 

Now look at the region where the argument of the theta function is positive. One 
gets two roots 



-p ± p^ + q^ - z"^ - 2ez - 

*i,2 = (.I'J 

z — q — e 



whose discriminant is negative for z ^ \—\Ip^ + 9^ + e, +£\/p'^ + + £]. Therefore, 

fe{q,p) = h{e\q,p) + h{e;q,p), (18) 

where 

2^3 f+oo p+00 1 
I,{e,q,p) = ^P^ dzze-^^'^ dt — , ^ ^ , r (19) 



t^'^Vt^ JVp^+e J -00 t'^{z-q) + 2pt + {z + q) 

and 



2a3 rV¥+¥+e ^-+00 j t2^,_ ,^^2pt + {z + q-e) 
h{e,q,p) = ^r^ dz dtO [ 



t^{z-q)+2pt + {z + , 



(20) 



Let us evaluate /i(e, q,p)- The integration over t yields ir/^/ z^ — q^ — p^ and thus 



2o? 

I,{e,q,p) = ^l , r , , dz. (21) 



7r\/7F J^/^i+^+e \/z^ - q^ - p2 
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An integration by part gives 

4c,5 r+oc 



= ^e-°'(9'+f') / yV°Vrfy + o(V^), (22) 

where = — — . Since the Gaussian integral equals y^/2a^, we finally get 

lim/i(e,g,p) = (23) 

eJ-O 

Therefore, it remains to prove that l2{£,q,p) vanishes for e — > 0. We will leave it 
as a very instructive exercise. 

Notice that the Radon transform can also be inverted by using the following 
alternative formula due to Helgason ([H), which is suitable for generalizations to 
symmetric homogeneous spaces 

/(g,p) = J-(-A)i/2 /■ fU^^qcos^ + psin^) d^. (24) 
4vr Jo 

Here the fractional Laplacian 

is defined by a Fourier transform 

{-A)'/'g{q,p) = [ {h + k^Y'^aiku k,)e^(^''^+^'^^)^Mh, (26) 



where 



g{ki,k2)= I g{q,p)e~^(<^''^^P'^^)^. (27) 



IR2 



In our case we would have to compute 

^ / q2 a2 \ 1/2 /•2tt 

" f _ _ \ / -a=^(gcosi?+psini3)2 j o 

ATT^y dq^ dp^J Jo 
a task even more difficult than the previous one. 



(28) 



4. Extension to n dimensions and discussion 



The definitions and conclusions of the previous sections can be easily extended to 
n dimensions. The Radon transform of a function / of the n-dimensional vector 



46 Journal of Modern Optics m2 final 



6 Taylor & Francis and I. T. Consultant 

x € IR" reads 



/»(u;,X)= / f{x)6{X -{oj,x))dx, 

JR" 



(29) 



where lij G S" ^, the (n — l)-dimensional sphere, (•, •) denotes scalar product and 
X £ R. The transform is 



f^'' {fi, X) = I f{x)5{X - (/., x)) dx, (30) 



AP 



where /i G IR" and X £ R. Obviously, from f {fj,,X) one can immediately recover 
P{uj,X) by setting fi = uj £ S""^: 

fHu^,X) = f''\u;,X). (31) 

However, notice that, although is the restriction of /^^" on the unit sphere 
S""^, there is actually a bijection between the two transforms. Therefore, they 
carry exactly the same information. Indeed, since the Dirac distribution is positive 
homogeneous of degree —1, i.e. 6{ax) = \a\~^6{x), for every a / 0, one gets from 
Eq. ^ 

fAP,.. v\ _ ^ fhP f 1^ ^ 



m VImI mJ 

for /U 7^ 0. In words, the tomogram /*'^^(/i, X) at a generic point ^ G IR" is com- 
pletely determined by the tomogram at G S""^. But the latter is nothing but 
the Radon transform, by Eq. (|3ip . Therefore we get the bijection 

/tt(^,X)=/*^^(^,X), (33) 

f''\^^^^) = hf'{^-y^-^ (^^0). (34) 

Notice also that at the origin fj. = the transform 

f''\0,X)=5{X) [ f{x)dx, (35) 

depends only on the total mass. 
The inversion formulae for the transforms (1291) and (1301) read 



2n^r. 

and 



/(^) = — 4rT(-A)("^^)/2 / /«(u;,(^,x)) du; (36) 



/*^^(/.,X)e^(^-(-))^, (37) 



respectively. While formula (j37p . which is nothing but a Fourier transform, is quite 
easy to handle, the inversion formula ()36p is in general very hard to tackle, espe- 
cially for even n, due to the presence of a fractional Laplacian. Therefore, from a 
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practical point of view our message is the following: in order to invert the Radon 
transform pip , dilate it by (I34p into the transform and then use the Fourier 
inversion formula ([37]). This simple trick enables one to avoid long and tedious 
calculations. 
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